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Abstract 

Non-supersymmetric SO(IO) grand unified theories provide a framework in which 
the stability of dark matter is explained while gauge coupling unification is realized. 
In this work, we systematically study this possibility by classifying weakly interact¬ 
ing dark matter candidates in terms of their quantum numbers of SU(2)i 0 U(l)y, 
B — L, and SU(2)/j. We consider both scalar and fermion candidates. We show that 
the requirement of a sufficiently high unification scale to ensure a proton lifetime 
compatible with experimental constraints plays a strong role in selecting viable can¬ 
didates. Among the scalar candidates originating from either a 16 or 144 of SO(IO), 
only SU(2)i singlets with zero hypercharge or doublets with P = 1/2 satisfy all con¬ 
straints for SU(4 )c«)SU(2)l«)SU(2)/j and SV{3)c^S\J{2)l0S\J{2)r0\J{1)b-l in¬ 
termediate scale gauge groups. Among fermion triplets with zero hypercharge, only 
a triplet in the 45 with intermediate group SU(4)c(8)SU(2)i(8)SU(2)j:j leads to solu¬ 
tions with Mqut > Afint and a long proton lifetime. We find three models with weak 
doublets and Y = 1/2 as dark matter candidates for the SU(4)c(8)SU(2)2,(8)SU(2)j:j 
and SU(4)c'(g) SU(2)2, (g) U(l)/j intermediate scale gauge groups assuming a minimal 
Higgs content. We also discuss how these models may be tested at accelerators and 
in dark matter detection experiments. 



1 Introduction 


Various cosmological observations have now established that more than 80% of the energy 
density of matter in the Universe is composed of non-baryonic dark matter (DM) [^. The 
Standard Model (SM) of particle physics, however, cannot explain this observation, and 
therefore there should be new physics beyond the SM which contains a DM candidate. 
One of the most promising class of candidates for DM is the so-called weakly-interacting 
massive particle (WIMP). These are electrically neutral and colorless particles which have 
masses of 0(10*^^“^^) GeV and couple to SM particles via weak-scale interactions. Their 
thermal relic abundance can explain the current energy density of DM. Such particles 
are predicted in many new-physics models; for example, the lightest neutralino in the 
supersymmetric (SUSY) SM is a well-known candidate for WIMP DM [^. 

For a WIMP to be DM, it should be stable or have a sufficiently long lifetime compared 
to the age of the Universe. To assure that, it is usually assumed that there is a symmetry 
which stabilizes the DM particle. For instance, in the minimal SUSY SM (MSSM), 
i?-parity makes the lightest SUSY particle stable and thus a candidate for DM in the 
Universe [^. Similarly, Kaluza-Klein parity in universal extra dimensional models 
and T-parity in the Littlest Higgs model yield stable particles, which can also be 
promising DM candidates. The ultraviolet (UV) origin of such a symmetry is, however, 
often obscure; thus it would be quite interesting if a theory which offers a DM candidate 
and simultaneously explains its stability can be realized as a UV completion rather than 
introducing the additional symmetry by hand. 

In fact, grand unihed theories (GUTs) can provide such a framework. Suppose that the 
ra nk of a GUT gauge group is larger than four. In this case, the GUT symmetry contains 
extra symmetries beyond the SM gauge symmetry. These extra symmetries should be 
spontaneously broken at a high-energy scale by a vacuum expectation value (VEV) of a 
Higgs field. Then, if we choose the proper representation for the Higgs held, there remains 
discrete symmetries, which can be used for DM stabilization [5 10 . The discrete charge 
of each representation is uniquely determined, and thus we can systematically identify 
possible DM candidates for each symmetry. 

In this work, we discuss the concrete realization of this scenario in non-SUSY SO(IO) 
GUT models. It is widely known that SO(IO) GUTs 


11 


13 have a lot of attractive 


features. Firstly, all of the SM quarks and leptons, as well as right-handed neutrinos, can 
be embedded into 16 representations of SO(IO). Secondly, the anomaly cancellation in 
the SM is naturally explained since SO(IO) is free from anomalies. Thirdly, one obtains 


improved gauge coupling unihcation |13| - |20| and improved fermion mass ratios 13 21 
if partial unihcation is achieved at an intermediate mass scale. In addition, since right- 
handed neutrinos have masses of the order of the intermediate scale, small neutrino masses 
can be explained via the seesaw mechanism if the intermediate scale is sufficiently 
high. SO(IO) includes an additional U(l) symmetry, which is assumed to be broken 
at the intermediate scale. If the Higgs held that breaks this additional U(l) symmetry 
belongs to a 126 dimensional representation, then a discrete Z 2 symmetry is preserved at 
low energies. One also hnds that as long as we focus on relatively small representations 
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(< 210 ), the 126 Higgs field leaving a Z 2 symmetry is the only possibility for a discrete 
symmetry 23,24 . We focus on this case in the following discussion. 


DM candidates appearing in such models can be classified into two types; one class of 
DM particles have effectively weak-scale interactions with the SM particles so that they 
are thermalized in the early universe, while the other class contains SM singlets which 
are never brought into thermal equilibrium. In the latter case, DM particles are pro¬ 
duced out of equilibrium via the thermal scattering involving heavy (intermediate scale) 
particle exchange processes. This type of DM is called Non-Equilibrium Thermal DM 
(NETDM) 1^, whose realization in SO(10) GUTs was thoroughly discussed in Ref. 24 


NETDM is necessarily fermionic as scalar DM would naturally couple to the SM Higgs 
bosons. Depending on the choice of the intermediate-scale gauge group, candidates for 
NETDM may originate from several different SO(IO) representations such as 45 , 54,126 
or 210 . Although the NETDM candidate itself does not affect the running of the gauge 
couplings from the weak scale to the intermediate scale, part of the original SO(IO) mul- 
tiplet has a mass at the intermediate scale and does affect the running up to the GUT 
scale. Demanding gauge coupling unification with a GUT scale above 10^® GeV leaves 
us with a limited set of potential NETDM candidates. When we further demand the 
splitting of the GUT scale or intermediate scale multiplets so that only a singlet survives 
at low energies only two candidates were left: in SU(4)c <S) SU(2)l 0 SU(2)ij notations 
there were the Dirac ( 1 , 1 , 3 ) originating in a 45 or a Weyl ( 15 , 1 , 1 ) also originating from 
a 45 or a 210 in SO(IO). 

In this paper, we study the remaining possibility, namely, that of WIMP DM can¬ 
didates in SO(IO) GUT models. We systematically classify WIMP DM candidates in 
terms of their quantum numbers and embed them in SO(IO) representations. As noted 
above and discussed in Ref. 24 , the presence of DM multiplets significantly affects the 


running of the gauge coupling constants. In this case, since the DM candidates are no 
longer SM singlets, the running of the gauge couplings is also affected between the weak 
and intermediate scales and may spoil gauge coupling unification realized in non-SUSY 
SO(IO) GUTs. We list WIMP DM models in which gauge coupling unification is achieved 
with appropriate GUT and intermediate scales. Then, we study the phenomenology of 
these models, such as the relic abundance of DM, the DM direct detection rate, the pro¬ 
ton decay lifetime, and neutrino masses. It is found that the condition of gauge coupling 
unification, as well as the proton decay bounds, severely restricts the WIMP DM models 
in SO(IO) GUTs. Still, we obtain some promising candidates, which can be probed in 
future DM searches and proton decay experiments. 

This paper is organized as follows. In the next section, we show the model setup for 
the SO(IO) WIMP DM scenario. The realization of a Z 2 symmetry and the classification 
of WIMP DM candidates are discussed there. Then, we analyze the scalar and fermionic 
DM models in Sec. and Sec. respectively. Section is devoted to conclusion and 
discussion. 
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2 Model 


We begin with an overview of the basic SO(IO) model needed to accommodate a DM 
candidate. As mentioned above, in this work, we consider SO(10) GUT models and 
restrict ourselves to a two step simultaneous symmetry breaking chainj^ in which the 
SO(IO) gauge group is broken to an intermediate gauge group Gmt at the GUT scale 
ATgut, and subsequently broken to the SM gauge group Gsm = SU(3)(^® SU(2)j;^®U(l)y 
and a Z 2 symmetry at the intermediate scale M^t- 


SO(IO) 


G 


int 


Gsm G) Z2 , 


( 1 ) 


The Higgs multiplets that break SO(IO) and Gjnt are labeled by Ri and R 2 , respectively. 
As discussed in the introduction, this Z 2 symmetry is a remnant of an extra U(l) symmetry 
in SO(IO) [5]-[^ and is used to stabilize DM candidates [9,10. A brief introduction to 

Possible SO(IO) 


2.1 


the intermediate subgroups and Z 2 symmetry will be given in Sec 
multiplets that contain an electric and color neutral component for a WIMP DM candidate 
are summarized in Sec. |2.2[ For a group theoretical argument on the classihcation of these 
DM candidates, see Appendix Among them, those who have a non-zero hypercharge 
are severely restricted by the DM direct search experiments. We consider this class of 
DM candidates in Sec. 12.31 and discuss conditions for the DM models to evade the direct 
search bound. 

To keep our model concise, in the following discussion, we only consider SO(IO) irre¬ 
ducible representations with dimensions up to 210. 


2.1 SO(IO) GUT and discrete symmetry 


We start by giving a brief description of the ingredients in our model. In an SO(IO) unifi¬ 
cation theory, a generation of SM fermions and a right-handed neutrino are embedded in a 
16 chiral representation, while the SM Higgs boson usually lies in a 10 representation. To 
obtain a realistic Yukawa sector, it is necessary to take the 10 to be complex 
will keep this sector unchanged in most of what follows. In addition to the SM particles, 
the Ri and R 2 Higgs representations are added to break SO(IO) and Gmt, respectively. 
The last ingredient of our model is the DM multiplet, whose lightest component is targeted 
to be the DM in the Universe. The stability of the DM is guaranteed by a remnant Z 2 
symmetry of the extra U(l) gauge symmetry of SO(IO) as we will discuss soon. Possible 
representations for the DM multiplet are determined below. Here, we assume that only a 
minimal set of the Higgs and DM multiplets which are necessary for the symmetry break¬ 
ing and mass generation of DM lie in the low-energy regime and other components have 
masses of the order of the symmetry breaking scale at which their masses are generated. 
For example, among the 10 representation, only the electroweak doublet components 
remains light to break the electroweak symmetry, while the other components have GUT- 
scale masses. Also, to obtain the right relic abundance, the mass of the DM particle is 


27,28 . We 


^For recent work on this kind of SO(IO) scenario, see Ref. 
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taken to be of order the TeV scale, while the masses of the other components are either 
of C>(Mgut) or C>(Mint). Such a hierarchical mass spectrum is obtained with hne-tunings 
similar to the doublet-triplet splitting needed for the 10 . In principle, it may be possible 
to achieve the splitting of the DM multiplets with a more elaborate scheme of particle 
representations as in the Higgs doublet-triplet separation (29}jM] . 

SO(IO) is a rank-hve group so it contains an additional U(l) symmetry besides the 
SM gauge symmetry. This additional U(l) can be broken into a Z 2 symmetry by a VEV 
of an appropriate Higgs held. If we restrict our attention to representation of dimension 
210 or less, the only choice of an irreducible R 2 that ensures the Z 2 symmetry is a 1260 
This Z 2 symmetry is equivalent to matter parity Pm = (— j^, under which the 
SM fermions are odd, while the SM Higgs held is even. Thus a fermion (boson) is stable 
if it has an even (odd) matter parity. 

We list in Table [T] all possible rank-hve subgroups and corresponding Ri whose VEV 
breaks SO(IO) to the subgroups. Here D denotes the so-called D-parity 33 , that is, a Z 2 
symmetry with respect to the exchange of SU( 2 ) 2 , O SU(2)j:j. D-parity can be related to 
an element of SO(IO) under which a fermion held transforms into its charge conjugate. 
In cases where D-parity is not broken by i?i, it is subsequently broken by R 2 at the 
intermediate scale M^t. In this work we only consider the subgroups without an explicit 
SU(5) factor. Since the DM is necessarily a color singlet, the running of the strong gauge 
coupling is unaltered by the presence of a new DM particle below the intermediate scale. 
Thus even though the addition of a DM multiplet yields unihcation of the gauge couplings, 
the unihcation scale Mint is always less than 10^“^ GeV as the contribution to the U(l)y 
beta function is always positive. If we now associate Mint with SU(5), this low partial 
unihcation is heavily disfavored on the basis of proton lifetime constraints. Flipped SU(5) 
usually has a high intermediate scale and a high GUT scale close to the Planck scale. In 
this case higher dimension operators suppressed by Planck scale become important, and 
one may also need to rely on a double seesaw for the explanation of neutrino masses |Mp5] . 
These bring extra complication into our model and we do not consider these possibilities 
here. Other intermediate gauge groups in the table are subgroups of SU(4)c' G SU( 2 ) 2 , ® 
SU(2)/j ® D, and U(l)^_^ is a subgroup of SU(4)c. The relationship among hypercharge 
V, the U(1 )b_£, charge B — L, and the third component of the SU(2 )k generators is 
very useful for determining the quantum numbers of DM candidates: 


Y 


B-L 


-R 


( 2 ) 


The convention we are using for hypercharge is such that electric charge is given by 
Q = Tl + Y, with T| denoting the third component of the SU(2)i;, generators. 

To summarize, our model contains the usual SM content and an SU(2)^ multiplet for 
DM at low energy scale. At the intermediate scale, depending on Gmt there are parts 
of the DM SO(IO) multiplet, parts of the 126 Higgs held to break while conserving 
matter parity, and perhaps some other Higgs helds that we specify on a model by model 


^The next-to-minimal possibility is a 1728. In addition, we note that a Z 3 symmetry can be obtained 
by a 672 Higgs field 
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Table 1: Candidates for the intermediate gauge group Gmt 

Gint 

Ri 

SU(4)c0SU(2)i0SU(2)^ 

210 

SU(4)c 0 SU(2)i 0 SU(2)fi 0 D 

54 

SU(4)c0SU(2)i0U(l)R 

45 

SU(3)^ 0 SU(2)i 0 SU(2)jj 0 U(l)B_i 

45 

SU(3)c 0 SU(2)i 0 SU(2)jj 0 U(1 )b_l 0 D 

210 

SU(3)c 0 SU(2)i 0 U(1 )r 0 U(l)s-L 

45 , 210 

SU(5) 0U(1) 

45 , 210 

Flipped SU(5) 0 U(l) 

45 , 210 


basis and are necessary for fine-tuning the DM mass. All other components in SO(IO) 
multiplets are assumed to be at the GUT scale. 


2.2 WIMP DM candidates 

In this section we discuss possible DM candidates in SO(IO) representations up to 210 , 
and classify them according to their quantum numbers. As discussed in the last section, 
the stability of DM is ensured by matter parity. Thus a fermionic DM candidate should 
be parity even and belong to a 10 , 45 , 54 , 120 , 126 , 210 or 210 ' representation, while 
scalar DM is parity odd and belongs to a 16 or 144 representation [8|^ 24 . Following the 
branching rules given in Ref. 36 , in Tablewe list SU(2)^ 0 U(l)y multiplets in various 
SO(IO) representations that contain an electrically neutral color singlet. A similar list of 
candidates can be found in earlier work . The table is classihed by R — L so one can 
check the matter parity of the candidates easily; R — L = 0, 2 candidates are fermionic 
while B — L = 1 candidates are scalar, labeled by an “F” or “S” at the beginning of each 
row, respectively. The subscript of the model names denotes the SU(2)i representation, 
while the superscript shows hypercharge. A hat is used for B — L = 2 candidates. 

We consider the WIMP DM scenario, which requires DM to be in thermal equilibrium 
with the SM particles before its abundance freezes out. This generally requires DM 
particles to interact with SM particles efficiently in the early universe. As a consequence, 
the fermionic singlets and F^ are not good WIMP candidates since they are SM singlets 
and can only interact with SM particles through exchange of intermediate scale virtual 
particles. In fact, these are examples of NETDM and the possibilities for F° and F° 
candidates for DM were discussed extensively in [^. Indeed, there it was shown that 
only two NETDM candidates from SO(IO) survive all phenomenological constraints. One 
possibility is associated with the SU(4)c 0 SU(2)i 0 SU(2)ij intermediate gauge group. 
In this case, the DM candidate is in a ( 1 , 1 , 3 ) originating in a 45 . This is an example 
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Table 2: List of SU(2)]^ ® U(1)y multiplets in SO(IO) representations that contain an 
electric neutral color singlet. 


Model 

B-L SU(2)^ 

Y 

SO(IO) representations 

F° 

1 

0 

45, 54, 210 

rp 

2 

1/2 

10, 120, 126, 210' 

pO 

3 

0 

0 

45, 54, 210 

pi 

3 

1 

54 

pl/2 

4 

1/2 

210' 


4 

3/2 

210' 

S? 

1 

0 

16, 144 

cl/2 

2 

1 

1/2 

16, 144 

cO 

^3 

3 

0 

144 

S3 

3 

1 

144 

F° 

1 

0 

126 

91/2 

2 2 

1/2 

210 

pi 

^3 

3 

1 

126 


of fJ. The second example is based on SU(4)c ® SU(2)i ® SU(2)/j 0 D and consists of a 
(15, 1, 1) originating from either a 45 or a 210 in SO(IO). Since the 15 of SU(4)c' carries 


lis is also an example of F^. All possible candidates associated with 
24 . A fermion that is a singlet under the intermediate gauge group 


zero B — L charge, t 
F^ were excluded in 

can also be produced through the exchange of the GUT scale particles, and thus be a 
DM candidate. For example, the case of the (1,1,1) component of a 210 is discussed in 
Ref. 


24 , which is again an example of F^ DM. 


The scalar singlet Sj and triplet Sg can interact with the SM Higgs boson efficiently 
through the quartic coupling and are potential good DM candidates to be discussed below. 
These can be taken to be either real or complex. For S^, there is no difference in any of 
our results whether Sj is real or complex. We have taken Sg to be real, but there would 
be no qualitative difference in our results for complex Sg. In addition, Sg couples to the 
SM particles via the weak interaction. Similarly, the fermion triplet Fg is a wino-like 
DM candidate and will also be considered below. In general, the neutral component of 
a SU(2)^ 0 U(l)y multiplet can interact with SM particles through exchange of kF or 0 
boson, and thus can be a good DM candidate. Such DM candidates have been widely 
studied in the literature 37-47 . 

F^l F 


There are also DM candidates which have non-zero hypercharge. These are: 12 


n/2 „3/2 


cl/2 
^2 1 


s^, F 


1/2 


3) 


and Fg. 


These DM candidates are severely constrained by DM 
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direct detection experiments since their scattering cross sections with a nucleon induced 
by Z-boson exchange are generally too large. Possible ways to evade this constraint are 
discussed in the following section. 


2.3 Hypercharged DM 


A DM candidate with V ^ 0 needs to be a Dirac spinor or a complex scalar, depending 
on its matter parity. These hypercharged candidates are severely restricted by the direct 
detection experiments, since they elastically scatter nucleons via the vector interactions 
mediated by Z-boson exchange, whose scattering cross section turns out to be too large by 
many orders of magnitude. One possible way to evade the constraint is to generate mass 
splitting. Am, between the neutral components of such a DM multiplet V’ and to split it 
into two Majorana fermions or real scalars xi, X 2 - Then, the DM no longer suffers from 
large scattering cross sections since it does not have vector interactions. Such splitting 
occurs if the DM mixes with extra SU(2)^®U(l)y multiplets after electroweak symmetry 
breaking, just like higgsinos in the MSSM, which originally form a Dirac fermion, reduce 
to neutralinos after they mix with gauginos. As we have assumed that only a single DM 
multiplet lies in the low-energy region, a natural mass scale for the extra SU(2)^ 0 U(l)y 
multiplets is the intermediate scale Mint. The effects of these heavy particles on the low- 
energy theory are expressed in terms of effective operators induced after integrating them 
out. Among them, the following operator is relevant for the generation of mass splitting 
for Dirac fermion DM: 




M-\ 

int 


(3) 


where H is the SM Higgs held, C represents charge conjugation, p = AY^ with > 0 being 
the hypercharge of the DM -0, and n = p — 1. For complex scalar DM, we have a similar 
operator with n = p — 2. Notice that the above operator violates any particle number 
assigned to the fermion 0. After the Higgs held acquires a VEV, the above operator 
reduces to a Majorana mass term, which generates a mass splitting between two Weyl 
fermions inside the neutral component of the DM multiplet. Namely, the mass eigenstates 
are expressed by two Majorana fermions, and the lighter one can be regarded as DM. The 
splitting is just given by Am ~ with v ~ 174 GeV being the Higgs VEV. Since a 

Majorana fermion does not couple to vector interactions, DM-nucleon elastic scattering 
cross sections are signihcantly reduced and we can avoid the direct detection bound. In 
the case of scalars, a similar operator to that in Eq. (|^ induces a splitting in the squared 
masses and hence Am ~ where m^ corresponds to the scalar mass term 

However, notice that the operator ([^ is considerably suppressed if the intermediate 
scale is large 42 . The suppression becomes more signihcant when IT is larger. In this 


case, the resultant mass splitting becomes extremely small. When the mass splitting is 
sufficiently small, then the DM can scatter off a nucleon N inelastically: XiA-N —)■ X 2 YN. 
Since this process is induced by the vector interactions, the scattering cross section again 
becomes too large if the mass splitting Am is smaller than the recoil energy. This sets 
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the bound Am > 100 keV. For this condition to be satisfied, ^ 10®, 3 x 10^, and 


4x 10^ GeV are required for fermionic dark DM with = 1/2, 1 and 3/2, respectively 42 
In the case of scalar DM, the upper bound depends on the DM mass. For a 1 TeV DM 
mass. Mint ^ 10® GeV for = 1. For a V;, = 1/2 scalar DM candidate, on the other 
hand, the mass splitting can be generated with a renormalizable interaction and its effect 
on the mass splitting depends only on its dimensionless coefficient. We will see later that 
this coefficient can still be very small, whose size is determined by the symmetry breaking 
pattern and its scale. This is because the operators relevant for the generation of the 
mass splitting are forbidden by the SO(IO) symmetry. Thus, the constraint from inelastic 
scattering can again give a bound on the intermediate scale even for = 1/2 scalar DM 
candidates. 

Another possibility to evade the direct detection bound is to push the DM mass 
sufficiently high. Since the local number density of DM is inversely proportional to the 
DM mass as the DM local energy density is fixed, the DM direct detection constraints 
are relaxed if the DM mass is taken to be heavy enough. In this case, DM is produced 
non-thermally 48 . This possibility is also discussed below. 


3 Scalar dark matter 


In this section, we discuss scalar WIMP DM in SO(IO) models with different intermediate 
subgroups. In this case, the DM candidates belong to either a 16 or a 144 representation. 
The masses of components in a DM multiplet in general need to be fine-tuned; if a charged 
component is nearly degenerate with the DM particle and decays to it only through an 
intermediate-scale gauge boson or Higgs field, this charged particle would be very long 
lived, which is cosmologically disastrous [^. Thus, to be safe, we take the masses of 
these extra components to be C>(Mgut) or C>(Mint), while the DM mass to be around 
TeV scale so that the thermal relic abundance of the DM agrees with the observed DM 
density, as we will see in Sec. 3.1 Here, we assume that the fine-tuning of DM masses be 


realized with a minimal choice of Higgs fields, that is, we exploit only Ri and = 126 to 
generate desired mass spectrum with Ri being an irreducible representation chosen from 
Table This is possible because a 16 or 144 can couple to the 126 Higgs field. Then, 
we study whether each set of matter content and its mass spectrum offers gauge coupling 
unification with appropriate GUT and intermediate scales. In Sec. 3.2, we present the 


results for the analysis and list promising candidates with M^t and Mgut determined 
by means of renormalization group equations (RGBs). The fine-tuning for the masses 
of components in a DM multiplet is discussed in Sec. |3.3[ As discussed in the previous 


section, hypercharged DM candidates require additional consideration for the generation 
of the mass splitting between the neutral components to avoid the bound from the direct 
detection experiments. This is discussed in Sec. |3.4[ Finally, in Sec. |3.5 , 


we summarize 


the current experimental constraints and future prospects for the scalar DM candidates 
discussed in this section. 











3.1 DM mass 


To determine the renormalization group (RG) running of gauge couplings, we need to 
know the mass of DM candidates, since they affect the running above its mass scale. An 
exception is as it is a SM singlet and does not contribute to the gauge coupling beta 
functions below Mint at the one-loop level. Scalar singlet DM is discussed in Ref. 49 . To 


roughly estimate favored mass region for such a singlet DM particle, consider the quartic 
interaction between the singlet DM (j) and the SM Higgs held: /2. Through 

this coupling, the singlet DM particles annihilate into a pair of the SM Higgs bosons. 
The annihilation cross section times the relative velocity between the initial state 
particles Urei is evaluated as 


^ann'i^rel — 


X2 




( 4 ) 


assuming that the DM mass moM is much larger than the SM Higgs mass rrih and we 
neglect terms proportional to The DM relic abundance is, on the other hand, related 
to the annihilation cross section by 




3 ^-1 


3 X 10 cm^ s 

(<^ann'f^rel) 


( 5 ) 


To account for the observed DM density = 0.12 [^, the DM mass should be 

^DM ^ 10 TeV for ^ 1- This gives us a rough upper bound for the DM mass. 

The other scalar DM candidates are SU(2)^ (g) U(l)y multiplets, which can interact 
with SM particles through gauge interactions besides the quartic coupling mentioned 
above. In particular, 83' is known as the Inert Higgs Doublet Model and has been widely 
studied in the literatur^ [^[^. To evaluate the effects of gauge interactions, let us hrst 
consider the limit of zero quartic couplings. In this case, the annihilation cross sections are 
completely determined as functions of the DM mass ruDM- Since the annihilation into SM 
fermions and Higgs boson suffers from p-wave suppression, the DM particles annihilate 
predominantly into the weak gauge bosons. In addition, we need to take into account 
coannihilation effects since all of the components in a SU(2)^ multiplet are degenerate in 
mass; the mass difference among the components is induced after electroweak symmetry 
breaking at the loop level and thus is quite suppressed compared to the DM mass, as 
small as 0(100) MeV. Taking these effects into account, for S^, the effective (averaged) 
annihilation cross section is given by 


C^'ann'i^rel — 


g^{3 - 4n2 + n^) + 16FV^ + 8g^g'^Y^{n‘^ 

647rc„m^j^ 



( 6 ) 


where g {g') are the SU(2)^ (U(l)y) gauge couplings and Cn = n (2n) for a real (complex) 
scalar. Here, we assume the DM mass to be much larger than the weak gauge boson 
masses. Again Eq. ([^ tells us that the masses of the DM candidates should fall into a 


^For another approach to the realization of the Inert Higgs doublet model based on grand unification, 
see Ref. 


52 
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region from ~ 500 GeV to ~ 2 TeV. On the other hand, if the quartic coupling is larger 
than the gauge couplings, the annihilation into a pair of Higgs bosons becomes dominant 
and thus the DM abundance would be similar to that of the singlet DM candidate. In 
general, the DM mass should lie between 0.5 TeV to 10 TeV for S2' , 


S3 and S3. 


More accurate estimations for the DM mass can be found in the literature 37,39 


53 with various additional contributions taken into account. For SU(2)^ ® U(l)y DM 
multiplets, the non-perturbative Sommerfeld enhancement is of great importance (54). In 


the limit of zero quartic coupling, the DM masses with which the relic abundance agrees 
with the observed DM density are evaluated as ttidm = 0.5 and 2.5 TeV for S2' and 
S3, respectively 53 . For S3, as far as we know, there has been no calculation which 


includes the Sommerfeld enhancement; thus its mass would be larger than 1.6 TeV, which 
is obtained with only the perturbative contribution considered 37 . For the cases where 


the scalar DM multiplets have non-zero quartic coupling with the SM Higgs doublet, it 
was shown in Ref. 39 that the allowed DM mass can be extended to ~ 58 and 28 TeV 
for S2' and S3, respectively, when the quartic coupling A ~ dvr. Such a large quartic 
coupling is, however, in general inconsistent with GUTs since it immediately blows up at 
a scale much below the GUT and intermediate scales. Thus, we implicitly assume the 
quartic coupling should be rather small, e.g., < 1, to avoid divergent couplings. In this 
case, the DM mass usually lies around 0{1) TeV. 


3.2 Candidates for scalar DM 

We list all possible scalar DM candidates in Table All of the candidates belong to either 
a 16 or a 144. Here, the hrst column shows the model names with subscript representing 
the intermediate gauge group Gint- The second column lists the Gint representations 
that contain the DM candidate multiplet S^. All of the components in the representation 
except the DM multiplet shown in the third column will have masses tuned to G(Mint). 
The rest of the components in the SO(IO) multiplet have masses of G(Mgut)- The case 
of Gint = SU(4)c (8 SU(2)i ® SU(2)fi ® D (SU(3)c ® SU(2)i ® SU(2)h 0 U(1)b_l 0 D) 
is identical to that of G^t = SU(4)c 0 SU(2)l 0 SU(2)r (SU(3)c 0 SU(2)l 0 SU(2)r 0 
U(l)B-i) with additional multiplets required by the left-right symmetry introduced above 
the intermediate scale. 

Gonsider, for example, SA422; the (4,1,2) DM multiplet originating from a 16 of 
SO(IO) in the SU(4)c’ 0 SU(2)i 0 SU(2)/j model. The “other half” of the 16, (4, 2,1), 
will have a GUT scale mass, while 7 of 8 helds in (4,1, 2) are tuned to have an intermediate 
scale mass and only the DM singlet is tuned to have a weak scale mass. In the case of 
SA422D, the DM multiplet corresponds to the (4,1, 2) © (4, 2,1) in the SU(4 )g 0 SU(2)i 0 
SU(2)^ 0 D model. In this case, none of the components have GUT scale masses and 
15 of the 16 fields have intermediate scale masses. Thus the presence of the left-right 
symmetry affects only the field content above the intermediate scale, though this will 
ultimately affect the scale of gauge coupling unification. These representations are added 
to a minimal SO(IO) unification model containing three generations of 16 chiral multiplet, 
a complex 10 Higgs multiplet, a. R2 = 126 Higgs multiplet, and a Ri Higgs multiplet 
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Table 3: Summary of DM multiplets. The second column shows the Gint representation 
with quantum numbers listed in the same order as the groups shown in the direct product. 
The case ofCi^t = SU{A)c®SUi2)L®SUi2)R®D (SUi3)c®SU{2)L®SU{2)R® U{1 )b-l® 
D) IS identical to that ofG^^t = SU{A)c® SU{2 )l® SU{2)r (SU{3)c® SU{2)l® SU{2)r® 
V(1)b-l) with additional multiplets required by left-right symmetry introduced above the 
intermediate scale. 


Model 

-Rdm 


SO(10) representation 

G 

int = SU(4)c ® SU(2)z. 

®SU(2)r(®D) 

SA422(D) 

4,1,2 

S° 

16, 144 

SB422(D) 

4,2,1 

ql /2 

16, 144 

SC422(D) 

4,2,3 

ql /2 

144 

SD422(D) 

4,3,2 

q1 

^3 

144 

SE422(D) 

4,3,2 

qO 

144 


Gint = SU(4)c 

®SU(2)i0U(l)H 

SA 421 

4,1,-1/2 

s? 

16, 144 

SB 421 

4,2,0 

ql /2 

^2 

16, 144 

SC 42 I 

4,2,1 

ql /2 

^2 

144 

SD 421 

4, 3,1/2 

q1 

^3 

144 

SE 421 

4, 3,-1/2 

qO 

144 

Gint = SU(3)c G SU(2) 

L 0 SU(2)ij (g) U(l)s L((g)D) 

SA3221(D) 

1 , 1 , 2,1 

S° 

16, 144 

SB3221(D) 

1 , 2 , 1,-1 

ql /2 

^2 

16, 144 

SC3221(D) 

1,2,3,-1 

ql /2 

^2 

144 

SD3221(D) 

1,3, 2,1 

q1 

^3 

144 

SE3221(D) 

1,3, 2,1 

qO 

^3 

144 


m 1 /2 m 

chosen from Table 1 Notice that 83 ' is contained in both the model SB’s and SC’s. The 
difference between the models is the SU(2)^ (or additional U(l)) charge assignment; for 
instance, SB422 (SC422) includes the SU(2)^ singlet (triplet) DM. From Table we hnd 
that a 16 contains only SA’s and SB’s, while a 144 has all of the candidates listed in the 
table. 

Next, we perform the RGEQ analysis in the models presented in Table |^to see if these 
models achieve gauge coupling unihcation with appropriate GUT and intermediate scales. 
The one-loop results for Mqut, ^int, the unihed gauge coupling acuT, and the proton 


^The beta functions for the minimal SO(IO) GUT described above are given in Appendix B of Ref. 


24 


11 

















Table 4: One-loop result for scales, unified couplings, and proton lifetimes for models in 
table. The DM mass is set to he moM =1 TeV. The mass scales are given in GeV 
and the proton lifetimes are in units of years. Blue shaded models evade the proton decay 
bound, t{p —)■ > 1.4 x 10^"^ 


yrs 


'55. 561 


Model 

logio ^GUT logic ^int 

CtGUT 

Tp{p e+TT^) 


Gint 

CO 

II 

®SU(2)iOSU(2)fl 

SA422 

16.33 

11.08 


0.0218 

36.8 ± 1.2 

SB422 

15.62 

12.38 


0.0228 

34.0 ± 1.2 

SC422 

14.89 

11.18 


0.0243 

31.0 ± 1.2 

SD422 

14.11 

13.29 


0.0253 

28.0 ± 1.2 

SE422 

14.73 

13.72 


0.0243 

30.4 ± 1.2 


Gint = 

SU(4)c ® SU(2)i 0 SU( 2 )r ® D 

SA422D 

15.23 

13.71 


0.0245 

32.4 ± 1.2 

SB422D 

15.01 

13.71 


0.0247 

31.6 ± 1.2 

SC422D 

14.50 

13.71 


0.0254 

29.5 ± 1.2 

SD422D 

13.95 

13.47 


0.0260 

27.3 ± 1.2 

SE422D 

14.55 

13.96 


0.0251 

29.7 ± 1.2 


Gint 

= SU( 4 )c 5 )SU( 2 )i®U(l)R 

SA421 

14.62 

10.96 


0.0226 

30.1 ± 1.2 

SB421 

14.55 

11.90 


0.0233 

29.8 ± 1.2 

SC421 

14.15 

10.92 


0.0236 

28.2 ± 1.2 

SD421 

13.91 

12.80 


0.0250 

27.2 ± 1.2 

SE421 

14.45 

13.12 


0.0241 

29.4 ± 1.2 


Gint = SU(3)c ® SU(2)i 

®SU(2)« 

® U( 1 )b l 

SA3221 

16.66 

8.54 


0.0217 

38.1 ± 1.2 

SB3221 

16.17 

9.80 


0.0223 

36.2 ± 1.2 

SC3221 

15.62 

9.14 


0.0230 

34.0 ± 1.2 

SD3221 

14.49 

12.07 


0.0246 

29.5 ± 1.2 

SE3221 

15.09 

12.22 


0.0237 

31.9 ± 1.2 


Gint = SU(3)c ® SU(2)i O SU(2)^ ® D{1)b-l ^ D 


SA3221D 

SB3221D 

15.58 

15.40 

10.08 

10.44 

0.0231 

0.0233 

33.8 ± 1.2 

33.1 ± 1.2 

SC3221D 

14.58 

11.62 

0.0245 

29.8 ± 1.2 

SD3221D 

14.07 

12.13 

0.0253 

27.8 ± 1.2 

SE3221D 

14.60 

12.29 

0.0245 

29.9 ± 1.2 
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lifetimes in the p —>■ e'^Ti^ channel are shown in Table [4|p| Here, Mqut and Mint are given 
in GeV units, while the unit for proton lifetimes Tp{p —)■ e~^7r^) is years. The DM mass 


is set to be tudm = 1 TeV. We have checked that altering the DM mass by an order 
of magnitude results in only a 0(0.2)% variation in the logarithmic masses of Mint and 
Mqut- The uncertainty of the lifetime reflects our innocence of the GUT-scale gauge 
boson mass Mx, which we take it to be within a range of O.SMqut ^ Mx < 2Mgut- 
It turns out that most models have already been ruled out by the current experimental 
constraint r(p —)■ > 1.4 x 10^^ yrs 55,56 . The models that possibly survive this 


constraint are SA422, SB422, SA3221, SB3221, SC3221, SA3221D, and SB3221D, which are highlighted 
in blue shading in the table. In terms of SU(2)2, G U(l)y assignments, only Sj and 83'^^ 
are found to be viable candidates. Among them, models SB422, SC3221, SA322 id 5 and SB3221D 
predict proton lifetimes close to the present limit, and thus can be tested in future proton 
decay experiments. 


3.3 Fine-tuning of scalar DM multiplets 

In the previous section, we have reduced the possibilities for Gmt to the only three gauge 
groups: SU(4 )cGSU(2)lGSU(2)k, SU(3)c0SU(2)i®SU(2)fi0U(l)B_L, and SU(3)cG 
SU(2 )L ® SU(2)ij G U(l) b-l ® D. According to Table = 210, 45, and 210 yield the 
above intermediate gauge groups, respectively. In this section, we briefly discuss how to 
obtain a desired mass spectrum for the DM multiplet using these i?i’s and i?2 = 126 with 
the help of hne-tuning. For convenience, we show an explicit procedure for the hne-tuning 
in Appendix by taking Rdm = 16 and Gjnt = 811(3)0 ® SU(2)o ® SU(2 )r 0 U(1)s_l 
as an example. 

Let ns hrst write down relevant terms for the mass terms of the DM multiplet A^dmI^ 


•^int 


= IVl l-ftDMl -r “T “r 

- + ^2\^Om\‘^\R2\‘^ + {Ai2® (1?Dm1?Dm)i26 (-^1" 

- Xi (1?dM-^Dm)45 {^1^1)45 + -^1 (-^DM-^Dm)210 (-^1-^1)210 

- (-RdM-^Dm )45 (-^21^2)45 + A2^° (-RdM-^Dm)210 (-^2-^2); 


210 ; 


(7) 


where the subscripts after the parentheses denote the 80(10) representation formed by 
the product in them. M, ki, and K 2 are dimensionful parameters, which we assume to 
be C>(Mgut)- Notice that the term (i?DM-RDM)i2o (-^1-^2)120 ^^*7 its charge conjugate 
vanish since the i?DM is a bosonic held and (AH) 120 is anti-symmetric with respect to 
the exchange of A and B. In addition, the term (i?DM7?DM)io (171-^2)10 Si^® ^ 

mass term for Hdm; {R 1 R 2 ) is singlet with respect to the 8M gauge interactions, and a 
10 representation does not contain such a component. The terms with the coefficients 

®We restrict our attention to one-loop running as two loop effects become very model dependent on 
our choice of the scalar potential. 

®In addition, there are couplings between the DM and the SM Higgs fields, which give a mass of the 
order of the electroweak scale to the DM multiplet. 
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A|; and are irrelevant to the generation of the mass splitting in the DM mnltiplet, as 
they only give a common mass to all of the components in the mnltiplet. It is also worth 
noting that terms inclnding K 2 and Aj|® break the particle nnmber which can be assigned 
to the complex scalar Rbm- Hence, these effects can split Rbm into two real scalars with 
different masses. We use these interactions to avoid the direct detection bound in the case 
of the complex hypercharged DM, models SB’s and SC’s, which we discuss in the following 
section. 

After Ri gets a VEV, the terms with ki, Af^, and A^^° generate mass terms for the 
components in Rbm with different mass values, since the Ri VEV couples to them with 
different Clebsch-Gordan coefficients. Thus, by hne tuning the coefficients M, ki, Af^, 
and A^^°, one can arrange that the DM mnltiplet obtains a mass of with other 

multiplets remaining around C>(Mgut)- 

The next step is to separate the 811(2)^;^ mnltiplet S^ from the intermediate gauge 
group mnltiplet. This can be accomplished by appropriately tuning the coefficients of K 2 , 
Aj;|®, Af^, and Af^° so that the generated mass terms cancel out the intermediate scale 
mass obtained previously, leaving only the DM candidate at TeV scaleAfter this step, 
we obtain a mass spectrum in which only the DM candidate lies around the TeV scale, 
while its partner helds with respect to the intermediate gauge symmetry are at M^t. The 
rest of the components of Rbm have masses of (9 (Mgut)- For an explicit example, see 
Appendix B. 


3.4 Mass splitting of hypercharged scalar dark matter 


As discussed in Sec. 2.3, we need a mass splitting of Am > lOOkeV 42 between the 


neutral and charged components of the hypercharged DM candidates (models SB and 
sc) to avoid the direct detection bound. Since both of these models yield S2' DM, the 
mass splitting can be induced by dimension-four operators like where (j) denotes 

the hypercharged scalar DM Such operators are, however, forbidden by the SO(IO) 
GUT symmetry; in fact, as the S^^^ DM and the SM Higgs held have B — L = 1 and 0, 
respectively, the operators contributing the mass splitting violate the B — L symmetry. 
Thus, they can be induced only below the intermediate scale where the B — L symmetry 
is spontaneously broken. 


Such an operator is induced by the interactions with the coefficients K 2 and Ajj 


126 


m 


Eq. ([^, since it requires violation of the particle number associated with the DM held 
(f). The required B — L breaking is realized by the R 2 VEV. We hnd that the tree-level 
diagrams in Fig. [T] are relevant to the generation of mass splitting. Here, Rh = 10 


^In the cases of models SB and SC, however, terms with K 2 and Aj;|® do not give mass terms for the 
DM mnltiplet after the intermediate symmetry breaking. The reason is as follows. For SB, since the i ?2 
VEV has a SU(2)ij charge while SB DM candidates do not, the couplings between the i ?2 VEV and the 
SB DM are forbidden by the SU(2)fl symmetry. For SC’s, since both the DM candidates and the i?2 VEV 
are in 811(2)^; triplets, a pair of the DM fields should be combined anti-symmetrically to be coupled to 
the i?2 VEV, but this vanishes because the DM is bosonic. In these cases, therefore, we carry out the 
fine-tuning for the DM mass only with the coefficients Af® and Af^^. 
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Figure 1: Diagrams that generate the mass splitting between real components of hyper- 
charged scalar DM. 


contains the SM Higgs field. Since the K 2 and interactions are symmetric with 
respect to the interchange of the component in R 2 which propagates in the upper 

two diagrams should be an 811(2)^ triplet. On the other hand, the component appearing 
in the inner lines in the lower two diagrams can be either an SU(2 )l singlet or triplet. 
The masses of these fields are dependent on the intermediate gauge groups; if G^t = 
SU(4)c 0 SU(2)i 0 SU(2)ij or SU(3)c:> 0 SU(2)i ® 811(2)/^ (g) U(1 )b_l, then these masses 
are C>(Mgut), while for Gjnt = SU(3)c ® 811(2)^ 0 SU(2)/j ® U(1)s_l ® -D, they are 
C>(Mint) because of the left-right parity D. 

Let us first consider the former cases. In these cases, the coefficient of the dimension- 
four operator is C>(Mint/MGUT), as the dimensionful couplings in the Lagrangian are 

expected to be (9 (Mgut)- We note here that there is no requirement for the cancellation 
between H 2 and to realize the desired mass spectrum since these couplings do 

not contribute to the mass splitting as mentioned in footnote in contrast to the case 
we discuss below. Thus, this operator induces a mass splitting of 


Am ~ 


^dmALgut 


( 8 ) 


The condition Am > 100 keV then becomes 


Mint 

Mgut 


> 3 X 10“® X 


/ ^DM \ 

VI Tevy 


(9) 


From Table we find that the model SB422 clearly satisfies this bound, while the inter¬ 
mediate scales in SB3221 and SC3221 he slightly below this constraint. However, since this 
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bound is just a rough estimation and the intermediate scales given in Table are obtained 
with the one-loop RGBs, it is possible that the DM candidates in these models are just 
not yet reached by the current direct detection experiments. If so, these models can be 
probed in the near future. 

For SB3221D, on the other hand, the mass spectrum is altered because of the pres¬ 
ence of the left-right parity. In this case, the charge of the DM candidate under 
is (1,2,1,—1), and the left-right symmetry requires the (1,1,2,-|-1) to be also light. 
To that end, the fine-tuning between the K 2 and Aj;|® terms in Eq. is required such 
that K 2 + ~ Mint; otherwise, these terms give a mass of O (a/M intM gut) to the 

(1,1,2,-|-1) component, which is much higher than the intermediate scale. This hne- 
tuning also guarantees the absence of non-perturbative couplings at low energies; without 
this fine-tuning, the exchange of intermediate-scale particles with the K 2 and Aj|®(i?i) 
vertices induces extremely large effective couplings, which destroy the perturbativity of 
the low-energy theory. 

In the presence of the fine-tuning, the diagrams in Fig. [^with the virtual states having 
a mass of Mint induce the effective operator with a coefficient of 0{1). Thus, the 

resultant mass splitting is well above 100 keV and the model SB3221D easily evades the 
constraints from the direct detection experiments. 

To summarize, SB422 and SB3221D are safe from the direct detection bound. SB3221 and 
SC3221 lie just around the margin of the bound, and they might be detected or completely 
excluded in future direct detection experiments. 


3.5 Constraints and prospects for the scalar DM candidates 

The above discussions have revealed that the only possible scalar DM candidates we could 
obtain with sufficiently high Mqut are S° and 83 ) as shown in Table [l Before concluding 
this section, we briefly review the current constraints on these DM candidates, and give 
prospects for probing them in future experiments. 

First, we consider Sj. This DM candidate has been widely discussed so far since it is 
one of the simplest extensions of the SM to include a DM candidate [^. As we have seen 
in Sec. 3.1, the thermal relic abundance of is determined once we fix the DM mass ut-dm 
and its quartic coupling to the SM Higgs field, Therefore, by requiring its thermal 
relic abundance to be equal to the observed DM density = 0.12, we can express 

\H(j> as a function of the DM mass tudm- Since this is the only coupling that connects the 
DM to the SM sector, various physical quantities relevant to the DM detection, such as 
the DM-nucleon scattering cross section, are also determined in terms of the DM mass. 

The present constraints on this DM are summarized in Refs. 57,58 . According to 


those results, currently, DM direct detection experiments give a stringent limit on the DM 
mass; the Sj DM with a mass of tudm ^ 53 GeV and 64 GeV < ruDM ^ 100 GeV has been 
excluded by the LUX experiment 59 . In addition, if tudm < with ~ 125 GeV 


the Higgs mass, then the constraint on the invisible decay width of the Higgs boson also 
restricts the DM. It turns out that the current upper bound on the invisible decay width 
BR(/;, —invisible) < 0.19 60 leads to the limit on the DM mass of moM ^ 53 GeV. The 
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DM direct detection experiments with ton-scale detectors, snch as XENONIT, will probe 



^ 100 GeV and moM ^ 500 GeV. In the former case, the DM particles annihilate 
efficiently throngh the Higgs boson exchange process, especially where ttidm — ''^h/2. 
When 100 GeV < ttidm ^ 500 GeV, the DM annihilation cross section is too large 
becanse the IV~^W~ channel is open. For ttidm ^ 500 GeV, both the Higgs boson and 
the gange bosons contribnte to the DM annihilation so that the resnltant relic abnndance 
can agree to the observed DM density. For the lower mass region, the direct detection 
experiments, the measurements of the Higgs decay branching ratios, and the electroweak 
precision measurements restrict the parameter space. Both of the mass regions can be 
probed in future direct detection experiments [^. Indirect detection experiments are, on 
the other hand, less promising; still, depending on the DM prohle, gamma-ray searches 
from the Galactic Genter may provide a signature of this DM candidate. 


4 Fermionic dark matter 


Next, we consider the fermionic DM candidates. 


Again, we begin with showing the 
As already mentioned above. 


favored mass region for these DM candidates in Sec. 4H 
the singlet fermion candidates, and F^, are not good candidates for WIMP DM since 
their annihilation cross sections are extremely suppressed (though they are good NETDM 
candidates). On the other hand, electroweakly charged DM can yield the desired relic 
abundance via gauge interactions. We discuss the V = 0 and V 7^ 0 cases in Sec. and 

Sec. 4^, respectively. We give some concrete examples for each case and perform RGE 
analysis to determine the intermediate/GUT scales of the models. Finally, in Sec 


4.4 


we 


summarize the present limits on these fermion DM models and discuss future prospects 


for probing these DM candidates. Non-thermal hypercharged DM is discussed in Sec. 4.5 


4.1 DM mass 

Gontrary to the case of the scalar DM, the thermal relic abundance of the fermionic DM 
candidates is completely determined by gauge interactions. Therefore, it is possible to 
make a robust prediction for the DM mass favored by the present DM density. In the 
case of fermion DM, not only the gauge boson channels but also the SM fermions and 
the Higgs boson hnal states can contribute to s-wave annihilation. We obtain a similar 
expression to Eq. (§ for the effective annihilation cross section of F^ as 

g^(2n^ + 17^2 - 19) + 4Y^g'^(41 + 8Y^) + 16^2^/2y2^^2 


<^ann'?^rel — 


128 irc„ 7 ii|,„ 


( 10 ) 


with c„ = 2n (4n) for a Majorana (Dirac) fermion. In addition, the Sommerfeld en¬ 
hancement again affects the annihilation cross section signihcantly. With this effect taken 
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into account, the thermal relic abundance of Fg is computed in Ref. 62 and found to 


be consistent with the observed DM density if tudm — 2.7 TeV as in the case of super 


37 


as m 
or the 


symmetric winos. As for and the favored mass value is ~ 1.1 TeV 
the case of supersymmetric Higgsinos. As far as we know, there is no calculation 
other fermionic DM candidates that includes the Sommerfeld enhancement; without the 
effect, the thermal relic of Fg, Fg, F^^, and is consistent with the observed value if 


ruDM — 1-9 TeV, 1.9 TeV, 2.4 TeV, and 2.6 TeV, respectively 53 


4.2 Real triplet DM 

We begin our discussion of fermionic DM models with the V = 0 case. As discussed earlier, 
these are less constrained by direct detection experiments. According to Table such 
candidates belong to SU(2)^ triplets in a 45, 54 or 210 of SO(IO). A summary of SU(4)c® 
SU(2 )l ® SU(2)r quantum numbers of these DM multiplets are listed in Table. Note 
that the B—L and charges for all of these DM candidates vanish, and therefore they are 
regarded as real Majorana fermions. As in the scalar DM scenario, the DM multiplet in the 
54 or 210 is degenerate with other components with respect to Gint, and we are required to 
break this degeneracy to avoid unwanted long-lived colored/charged particles 24 . In the 


fermionic case, however, a renormalizable Yukawa term like i?DMRDMl26// is forbidden 
by SO(IO) symmetry and the choice of DM representation 24 , and thus we are unable 


to use the 126 Higgs to break the degeneracy. Therefore, we need to introduce additional 
Higgs fields at the intermediate scale in these cases. 


Table 5: Real triplet DM eandidates in various SO(10) representations. 


SO(10) representation 

SU(4)^® SU(2)^0SU(2)^ 

45 

(1,3,1) 

54 

(1,3,3) 

210 

(15,3,1) 


For simplicity, we restrict ourselves to the cases where the intermediate scale VEVs 
develop in the SM singlet direction of Ri and/or R 2 = 126. One of the SM singlet 
components of Ri should have a VEV of C>(Mgut) to break SO(IO) into Gint. The 
R 2 Higgs held acquires an VEV to break Gjnt, but it is not able to give mass 

differences among the components in Rbm, as mentioned above. Thus, we need to exploit 
an extra SM singlet component in Ri which remains light compared to the GUT scale, 
to induce intermediate-scale mass terms for Rbm, which are to be used to generate the 
required mass splitting. We denote the VEVs of these two components of Ri which break 
SO(IO) and Gmt by ugut ~ Mqvt and Uint ~ Mint, respectively. Then, the mass splitting 
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in the DM multiplet -Rdm can be realized as follows: 

= d^-RoM-RoM — -Ri-RdM-RdM 

(M - CiWgut - C2ni„t) XX , (H) 

where x denotes the DM held and M ~ Mqut is a universal mass. Ci and C2 are the 
Clebsch-Gordan coefficients that vary for different Rdm components. Thus, by hne-tuning 
M such that M — ciUgut — C 2 V\nt ~ 1 TeV, we can set the DM triplet to be at TeV scale 
while leaving other contents in Rbm either around Mint or Mgut- We summarize in Table 
the multiplets in Ri that may develop a VEV of (9(Mint) for different Gint- The multiplets 
are labeled by the quantum numbers of Gint. It turns out that there is no extra SM singlet 
component in 54, which is indicated by a hyphen in the table. As a consequence, there 
is no way to hne-tune the mass of the (1, 3, 3) DM candidate originating from the 54 
and we drop it from further discussion. Here, we note that the cases of Gint = SU(3)^ G) 
SU(2 )l®SU(2)h(8)U(1)b_l and SU(3)f;0SU(2)L(8)SU(2)i?®U(l)B_L®D are disfavored 
before further analysis: the addition of a real triplet DM lowers the unihcation scale to 
unacceptable values and in these cases there is neither any new-physics contribution to 
the 811(3)0 gauge coupling beta function nor any new positive contribution to the 811(2)^ 
beta function above Mint. Therefore, the resultant Mgut is always smaller than the 
unihcation scale of the 8U(3)c and 8U(2 )l gauge couplings in the 8M plus a real triplet 
DM, which is below 10^® GeV and thus too low to evade the proton decay constraint]^ 
For this reason, we do not consider these cases in Table 


Table 6: Possible components in Ri that can develop a VEV of 0{Mint)- 


Gint 

Ri 

Intermediate scale multiplets 

8U(4)c®8U(2)d®8U(2)k 

210 

(15,1,1) 



(15,1,3) 

8U(4)c (8) 8U(2 )d ® 8U(2)k (g) D 

54 

- 

8U(4)c®8U(2)d®U(1)k 

45 

(15,1,0) 


We now perform the RG analysis to look for promising models with additional inter¬ 
mediate Higgs multiplets given in Table The one-loop results for Mgut, Afint, acuT, 
and proton lifetimes for different combination of Rdm and the Higgs helds are listed in 
Table [7jp| Here, we set the DM mass to be 1 TeV. The second column lists the extra Higgs 
helds in Ri at M^t in addition to R 2 . We suppressed combinations of Higgs multiplets 

®Note that scalar doublet DM is allowed under these intermediate symmetries as shown in Table 
since its contribution to the beta functions is much smaller than that from a fermionic real triplet, thus 
allowing for a higher unification scale. 

®We again restrict our attention to one-loop RGBs to avoid any model dependence due to the Yukawa 
coupling with the additional Higgs in Ri. 
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Table 7: The one-loop results for Mcut, Mint, chgut, and proton lifetimes for real triplet 
fermionic DM models. Here we set the DM mass to he 1 TeV. The mass scales and proton 
decay lifetime are in unit of GeV and years, respectively. In the blue shaded model, gauge 
coupling unification is achieved with a sufficiently high GUT scale. 


Rdm 

Additional Higgs 

in Ri 

logioMint logioMcuT 

ttGUT logiorp(p-)■ e+7r°) 


Gint 

= SU(4)^ 

0 SU(2)^® SU(2)^ 

(1,3,1) 

- 

15.50 

13.69 

0.0263 

(1,3,1) 

(15,1,3) 

- 

- 

- - 

(1,3,1) 

(15,1,1) 

15.65 

13.47 

0.0263 

(1,3,1) 

(15,1,1) 

(15,1,3) 

6.54 

17.17 

0.0252 39.8 ±1.2 

(15,3,1) 

(15,1,1) 

14.44 

14.10 

0.0246 

(15,3,1) 

(15,1,1) 

(15,1,3) 

14.52 

14.11 

0.0243 


Gint = 

SU(4)c ® SU(2)^ ® SU(2) 

r^D 

(1,3,1) 

- 

14.78 

14.04 

0.0250 


Gint 

= SU(4)c 

SU(2)^®U(1)^ 

(15,3,0) 

(15,1,0) 

14.55 

14.21 

0.0246 


that cannot split the degeneracy of DM multiplet as in Eq. (0. The mass scales and 
proton decay lifetime are in nnits of GeV and years, respectively. We hnd that there is 
only one promising model with Gjnt = SU(4)^ ® SU(2)^ ® SU(2)^, which is highlighted 
by blue shading in Table. In this case, since the DM multiplet is a singlet under both 
SU(4)(^ and SU(2)^, the additional Higgs helds are not necessary from the viewpoint of 
the mass splitting for the DM multiplet; namely, there is no degeneracy problem for this 
model. Rather, they are required so that the model achieves a good unihcation scale 
beyond proton decay constraint. The model has, however, a quite low intermediate scale 
that results in large neutrino masses through the type-I seesaw mechanism since the Dirac 
mass terms for neutrinos are related to the up-type Yukawa couplings in this setup. A 
simple way to evade this problem is to introduce a complex (15, 2, 2)c Higgs held in a 126 
to modify the relation, as discussed in Ref. 24 |Hlfa(15,2,2)c Higgs is also present at 


the intermediate scale, it turns out that gauge coupling unihcation is still realized, with 
logioMnt = 9.28, logioMcuT = 16.38, ogut = 0.038, and \og^QTp{p -)■ e+7r°) = 35.9. 
Here again, the mass scales and proton decay lifetime are expressed in units of GeV and 


10 


For the effects of a (15, 2, 2)c Higgs field on the Yukawa couplings, see Refs. 27 63 . 
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years, respectively. Finally, we note that the addition of (15, 2, 2)c will not resurrect the 
failed models in Table 0 


4.3 Hypercharged DM 


Hypercharged DM is a natural step forward after considering real triplet DM. In this 
section, we still restrict the Higgs content as in the previous section. As we discussed in 
Sec. |2.3[ hypercharged DM is strongly constrained by direct detection experiments. To 


evade this constraint, we need to split the mass of the Weyl components of the hyper¬ 
charged Dirac DM by ~ 100 keV. There are two possible ways to generate an effective 
operator in Eq. ([^ through exchange of a field at the intermediate scale at tree level, de¬ 
pending on whether it is a scalar or a fermion. In the former case, the effective operator is 
induced by the exchange of intermediate-scale Higgs fields, as illustrated in Fig. 2(a) This 


requires the hypercharge of the virtual Higgs held to be at least one and Mint ^ 10^ GeV. 
According to Table the only candidate for such a Higgs held belongs to (15,1,3) in 
the 210 when Gmt = SU(4)c 0 SU(2)i 0 811(2)^. The DM candidate should then be in 
a (15,2,2) or (10,2,2) © (10,2,2) representation of SU(4)c © SU(2)2, © SU(2)7^. We 
performed a scan for models that contain above content, and found that none of them 
gives appropriate Mint and Mqut. The latter possibility is to introduce another fermionic 
real multiplet at the intermediate scale, so that the DM candidate is a mixture of a hyper¬ 
charged held and a Majorana held. This mechanism is demonstrated in Fig. 2(b)| where 
Rbm is the main component of the DM candidate which is hypercharged and has a mass 
term of TeV scale; i^ the Majorana held at the intermediate scale. The cross mark 


in Fig. 2(b) represents the chiral hipping in the propagator of the Majorana held R[ 


Rbm and couple to the SM Higgs held through terms like 

.^mix oc RbmR'bm^h + h.c. 


DM- 


( 12 ) 


Since is a Majorana held, it can only belong to either a singlet or a real triplet 
among the possible candidates in Table As a result, DM can only belong to a doublet 


or or a quartet (F^^^), w ith hypercharge 1/2. This requires M^t < 10® GeV 




n/2 


according to the discussion in Sec. 


2.3 


Note that the V > 1 DM candidates, Fg, Fg, 


3/2 

and F4 , require at least 2Y additional fermions at the intermediate scale to generate the 
ehective operator in Eq. ([^. To minimize our model content, we do not consider these 
possibilities in the following discussion. 

Taking the above discussion into account, we list the possible SO(lO) representations 
for Rbm in the upper part of Table the singlet and real triplet candidates for are 
listed in the lower part of Table and Table respectively. The quantum numbers of the 
DM candidates with respect to the intermediate gauge groups we consider can be inferred 
from the SU(4)^ © SU(2)^ © SU(2)^ and B — L quantum numbers listed in the table. 

Now, we perform a one-loop calculation of Mint, Mqut and the proton decay lifetime 
for various combination of i?DM, -Rdm intermediate scale Higgs helds. Then, we pick 
up the models that are not yet ruled out by proton decay experiments, and at the same 
time have a relatively low intermediate scale M^t < 10®. We also require that the models 
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(a) Scalar exchange (b) Fermion exchange 

Figure 2: Diagrams that generate the mass splitting between the Weyl components of hy- 
percharged Dirac DM through the exchange of an intermediate-scale (a) scalar (b) fermion. 


Table 8: The upper half of the table shows the fermionic F = 1/2 candidates for Rdm 
in various SO(IO) representations; the lower half of the table shows the fermionic singlet 
candidates for 


SO(IO) representation 

SU(4)^®SU(2)^0SU(2)^ 

B-L 

10, 120, 210' 

(1.2,2) 

0 

120, 126 

(15,2,2) 

0 

210 

(10, 2, 2) ©(TO, 2, 2) 

±2 

210' 

(1.4,4) 

0 

54, 210 

(1,1.1) 

0 

45 

(1,1,3) 

0 

45, 210 

(15,1,1) 

0 

210 

(15,1,3) 

0 

126 

(10,1,3) 

2 


have appropriate particle and Higgs content, so that the DM acquires the right mass 
through Eq. (11) and Eq. (12). It turns out that the viable models are limited to Gjnt = 
SU(4)(^®SU(2 )j;^( 8)SU(2)^ or SU(4)(^®SU(2)^®U(l)j:j. These models are listed in Table|^ 
and no quartic models (F^^) were found. The model FA422 is incompatible with small 
neutrino masses, since the Yukawa coupling for the 16 of this model is unihed at Mgut- 
For models FA421 and FB422 , on the other hand, we can avoid the neutrino mass problem by 
hne-tuning the Yukawa couplings with additional Higgs helds at the intermediate scale, as 
discussed in Sec. 4Y Among them, the model FA421 has a phenomenologically interesting 
consequence. Since Mint — 3 TeV, this model predicts a new massive neutral gauge 
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boson, Z\ and vector leptoqnarks whose masses are aronnd a few TeV. These particles 
can be probed in fntnre LHC experiments; for instance, dilepton resonance searches 64 


are powerfnl probes for snch a Z'. The leptoqnarks are pair prodnced at the LHC, and 
their signatnre is observed in dijet plus dilepton channels 65 . Since they are produced 


via the strong interaction, their production cross section is quite large. Thanks to the 
distinct hnal states and large production cross section, the LHC experiments can probe 
TeV-scale leptoqnarks at the next stage of the LHC running. 


Table 9: Possible hypercharged fermionic DM models that is not yet excluded by current 
proton decay experiments. The quantum numbers are labeled in the same order as Gmt- 
The subscripts D and W refer to Dirac and Weyl respectively. The numerical results are 
calculated for DM mass of 1 TeV. The mass scales and proton decay lifetime are in unit 
of GeV and years, respectively. 


Model 

Adm 

Adm 

Higgs logioMnt logioMcux 

ogut 

logio A 



Dint 

= SU(4)ci® SU(2)^®U(1)^ 



FA421 

(l,2,l/2)t3 

(15,l,0)tt^ 

(15,l,0)i{ 3.48 17.54 

(15,2,1/2)c 

0.0320 

40.9 ± 1.2 



Gint^ 

= SU(4)^ (g SU(2)^ (g SU(2)^ 



FA422 

(1,2,2 V 

(l,3,l)w 

(15,l,l)fl 9.00 15.68 

(15,l,3)fl 

0.0258 

34.0 ± 1.2 

FB422 

(1,2,2 V 

(1,3,l)w. 

(15.1.1) fl 5.84 17.01 

(15.2.2) c 

(15.1.3) fl 

0.0587 

38.0 ± 1.2 


To conclude this section, we perform a scan for more general models where the addi¬ 
tional intermediate scale Higgs helds are not restricted to the ones in Ri. Instead, they can 
be any combination of Gjnt representations that contain SM singlets. The Higgs helds can 
be taken to be either real or complex. Moreover, we also consider the possible addition of a 
(15,2,2)c Higgs at the intermediate scale, which can be used to evade the problem of large 
neutrino masses. The result of the scan is demonstrated in a scatter plot in Fig. The DM 
mass is again hxed to be 1 TeV. The real triplet DM, doublet-singlet mixing 

DM and doublet-triplet mixing DM cases are colored in red, blue and green, respectively. 
The triangle, circle and square marker corresponds to G-^t = SU(4)^ ® SU(2)^ ® U(l)^, 
SU(4)^ 0 SU(2)^ O SU(2)^ and SU(4)^ O SU(2)^ O SU(2)^ 0 D, respectively. The 
Mint > ^GUT region is theoretically disfavored, and is indicated by the gray shaded 
area. In this plot, we do not consider the realizability of the mass hierarchy for the DM 
multiplet, and thus the number of good models should be smaller than that shown in the 
plot. All of the SU(4)(^ (g) SU(2)^ (g) SU(2)^ (g) D cases with doublet DM predict the same 
Mint, since the addition of extra helds at the intermediate scale does not change Mint 
in the presence of the left-right symmetry [^. As can be seen from the hgure, model 
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Figure 3: A scatter plot of general fermionic DM models. Here, the DM mass is set 
to he 1 TeV. The real triplet DM, doublet-singlet mixing DM and doublet-triplet mixing 
DM cases are eolored in red, blue and green, respectively. The triangle, cirele and sguare 
marker corresponds to Gint = SU(4)(^ 0 SU(2)^ ® 811(4)^ ® SU(2)^ ® SU(2)^ 

and SU(4)p 0 SU(2)^ 0 SU(2)^ 0 D, respectively. The vertieal dashed line at 10® GeV 
indicates the direct detection constraint for F = 1/2 dark matter. The gray shaded area 
is disfavored for M^t > Mgut- Only the filled symbols are consistent with a sufficiently 
long proton lifetime. 


points are concentrated in the high intermediate scale and low GUT scale region. After 
we apply the constraint of proton decay lifetime (shown by the hlled symbols) as well as 
the condition Mint ^ 10® GeV for the doublet DM cases, the viable candidates turn out 
to be fairly limited. 


4.4 Constraints and prospects for the fermion DM candidates 


Finally, we review the present experimental constraints and future prospects of the fermionic 
DM candidates discussed in this section. Let us begin with the real triplet DM case. At 
the LHC, this DM candidate can be probed by searching for disappearing tracks caused 
by the charged component of the triplet DM, which has a decay length of 0{1) cm. 
Such a small decay length is due to the small mass difference between the neutral and 
charged components; it is as small as a hundred MeV, since it is induced at loop levep] 


Currently, the mass difference is computed at two-loop level 
difference is about 165 MeV. 


66 : for a 3-TeV triplet DM, the mass 
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after electroweak symmetry breaking. Based on this search strategy, the ATLAS ex¬ 
periment has searched for triplet DM and has given a lower bonnd on the DM mass of 

For fntnre prospects on collider searches of triplet DM, see Ref. 1^. 


iriDM > 270 GeV 67 


Indirect searches of triplet DM are also promising since this DM has a large annihilation 
cross section, as already seen in Sec. 4.1 Indeed, an excess of cosmic-ray antiprotons 
observed by the AMS-02 experiment (Ml might be the first signatnre of triplet DM 70 


On the other hand, this DM is cnrrently being constrained by the searches for gamma-ray 
line spectrnm coming from the Galactic Genter. As discussed in Ref. [^, the results from 
the H.E.S.S. Gollaboration 72 may give a strong limit on triplet DM, though the conse¬ 


quences depend on the DM density prohle used in the analysis. More robust constraints 
are obtained by means of the observation of gamma rays from dwarf spheroidal galaxies 
given by the Fermi-LAT collaboration 73 ; according to that result, the mass of triplet 


DM is limited to 320 GeV < tudm < 2.25 TeV and 2.43 TeV < tudm at 95% confidence 
level 74 . In the future, gamma-ray search experiments can probe a wider range of masses 


region for triplet DM. Direct detection experiments are also able to catch the signature of 
triplet DM. The scattering of triplet DM with a nucleon is induced by the exchange of the 
electroweak gauge bosons at loop level 75 , and its scattering cross section is evaluated 


at the next-to-leading order in Ref. 76 


for instance, cxsi ~ 2 x 10 cm"^ for a 3 TeV 


real triplet DM, which is well above the neutrino background 
see also Ref. 


77 . For relevant works. 


78 . As a consequence, the triplet DM scenario can be tested in various 


future experiments, and therefore is a quite interesting possibility among the SO(IO) DM 
candidates. 

Next, we consider the doublet DM case. Gontrary to triplet DM, doublet DM is rather 
hard to probe in experiments. In this case, the mass difference between the neutral and 
charged components is found to be as large as several hundreds of MeV, which makes it 
difficult to search for doublet DM signal events at the LHG by using the disappearing track 
method. The most promising way to probe doublet DM is the direct pair production at the 
ILG, which also enables ns to study its properties precisely 79 . The indirect DM searches 


are also less promising due to a relatively small annihilation cross section of doublet DM. 
The direct detection of this DM is only possible when the intermediate scale is as low 


as 10^ GeV 42 ; in this case, the DM-nucleon scattering occurs through the exchange of 


the Higgs boson, which is induced by effective operators generated at the intermediate 
scale, ^ In addition, if there are additional GP phases in the effective operators, the 

For 


41 42 


electric dipole moment of electron is sensitive to the effects of doublet DM 
instance, the model FA421 and FB422 in Table may be tested in these experiments. After 
all, the prospects for probing doublet DM quite depend on the intermediate scale, and 
future experiments are able to search for some of the doublet DM models discussed above. 


^The electroweak loop contribution to the DM-nucleon scattering in this case turns out to be very 


small 76 . 
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4.5 Non-thermal Hypercharged DM 

In our previous discussion of fermionic hypercharged DM, we have assumed the DM to be 
in thermal equilibrium before freeze out so that the DM mass is restricted to be 0(1 TeV) 
by the observed abundance. As a consequence, hypercharged DM is highly restricted by 
direct detection experiments. One way around such a constraint is to introduce a small 
mass splitting between the Weyl components of the DM by mixing the DM with another 
Majorana multiplet at the intermediate scale. On the other hand, as suggested in Ref. (4^ , 
non-thermally produced DM can be extremely heavy and thus avoid the direct detection 
constraint even when there is no mass splitting between the Weyl components. As a 
result, a minimal hypercharged DM model is possible in this scenario. This minimality 
motivates us to consider this class of DM candidates, even though they are not WIMPs. 

First, let us derive the lower bound on the DM mass to evade the direct detection 
constraints. A fermionic DM particle with non-zero hypercharge Y scatters off a nucleus 
via the exchange of Z boson, and its scattering cross section is given by 


GlY'^ 


a = 


2ti 


[N - (1 - A sin^ 9w)Zf 


(- 


m-DM^T 


\mBM + rriT 


( 13 ) 


where Gj? is the Fermi constant, 6iy is the weak mixing angle, Z and N are the numbers 
of protons and neutrons in the nucleus, respectively, and rriT denotes the mass of the 
target nucleus. The LUX limit then readt 


13 


mBM > 


(2Yf X 6 X 10^ GeV . 


(14) 


Such a heavy DM candidate can lead to the correct relic abundance if its mass is larger 
than the reheating temperature Tr after inflation so that it never reaches equilibrium with 
the thermal bath of SM particles. Then by carefully choosing the DM mass, the reheating 
temperature, and the maximum temperature after inflation, one obtains the desired relic 
abundance. The reheating temperature was shown to be in a range of 


48 


Tr ~ 


moM 


3 X IQio GeV 


( 15 ) 


In this scenario, it is natural to assume the DM mass scale to be equal to intermediate 
scales of the unihcation models that we consider. Usually, Mint is large enough to evade 
the direct detection bound and yet not so large that gravitational production of DM 
becomes dominant. Moreover, we do not need to worry about mass splitting in the DM 
multiplet, as we did in the WIMP scenario. For such heavy particles, the electroweak and 
strong corrections to the masses of charged or colored particles in the multiplet are large 
enough to prevent them from acquiring cosmological lifetimes. This allows us to consider 
minimal hypercharged DM models, which contain only one DM multiplet without any 
extra Higgs multiplets with respect to the minimal SO(IO) unihcation model. 

lapoj. 131^6 target, Z = bA,N = 77. 
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The y = 1/2 DM candidates we consider inclnde those in the npper part of Table 
and the Y = 1 candidates listed in Table IT The last mnltiplet in the npper part of Table^ 
also contains an D = 3/2 candidate of DM which we will also consider. Now that the DM 
mnltiplet only contribntes to the rnnning above Mint, we will also consider the cases of 
(Tint = SU(3)^®SU(2)i®SU(2)^®U(l)B_i and SU(3)^®SU(2)i0SU(2)H®U(l)s_L®D, 
in contrast to the WIMP scenario. 


Table 10: Fermionic Y = 1 candidates for Rum in various SO(IO) representations; 


SO(IO) representation 

SU(4)^®SU(2)^®SU(2)^ 

B-L 

54 

(1,3,3) 

0 

126 

(10,3,1) 

2 


After an exhanstive calcnlation for different choices of Rbm and Gint, we hnd several 
possible models listed in Table [T^ that snrvive the direct detection and proton decay 
constraints. Viable minimal models only exist for donblet DM when Gjnt = SU(4)(^ ® 
SU(2)^ 0 SU(2)^, (Tint = SU(4)^ ® SU(2)^ 0 SU(2)^ 0 D and Gint = SU(3)^ 0 SU(2)^ ® 
SU(2)j:j (g) U(l)^_^. The intermediate scale of each model is larger than 10® GeV, large 
enongh to evade the direct detection bonnd indicated by Eq. (14). As can be seen, most 
of the models in Table IT can be probed in fntnre proton decay experiments, though it is 
hard to detect these DM candidates in direct detection experiments. 


Table 11: Possible non-thermal hypercharged fermionic DM models that are not yet 
excluded by current direct detection and proton decay experiments. The quantum numbers 
are labeled in the same order as Gint- The numerical results are calculated for DM mass 
equal to Mint- The mass scales and proton decay lifetimes are in unit of GeV and years, 
respectively. 


Model 

Rom 

logtoMnt 

logigMcux 

OguT 

logic A 


Gint 

= SU(4)c. ® SU(2)^ SU(2)^ 


FNA422 

(1,2,2)h. 

12.10 

15.63 

0.0225 

34.0 ± 1.2 

FNB422 

(15,2,2)h. 

11.15 

16.77 

0.0387 

37.9 ± 1.2 


Gint = 

SU(4)p® SU(2)^® SU(2) 



FNA422D 

(15,2,2)h. 

13.71 

15.36 

0.0286 

32.8 ± 1.2 

FNB422D 

(10,2,2)t, 

13.71 

15.94 

0.0342 

34.9 ± 1.2 


Gint = SU(3)^ ® SU(2)^ 

®SU(2)^®U(1)5_^ 


FNA3221 

(1,2,2,0)vk 10.34 

15.82 

0.0227 

34.8 ± 1.2 
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5 Conclusion and discussion 


The success of the Standard Model is now well established. Nevertheless, we know that 
the Standard Model is incomplete. Neutrinos have masses, there is a non-zero baryon 
asymmetry in the Universe, and dark matter makes up a sizable component of the total 
matter density. Many extensions to the Standard Model have been studied to explain 
these phenomena. But rarely can a single extension explain all three. SO(IO) grand 
unification is one such example. 

In most models, SO(IO) symmetry breaking involves an intermediate gauge group, 
whose unknown scale presumably lies between the weak scale and the grand unihed scale 
(defined by the renormalization scale where the gauge couplings are equal). Standard 
Model fermions are neatly contained in a 16 of SO(IO) which includes all of the known 
SM fermions plus a right-handed neutrino per generation. As the right-handed neutrino 
is a SM singlet, it easily picks up a mass of order the intermediate scale during the second 
phase of symmetry breaking, and leads to a natural realization of the see-saw mecha- 
for the generation of neutrino masses. If produced (thermally or non-thermally) 


nism 
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after inflation, these same right-handed neutrinos can decay out of equilibrium and pro¬ 
duce a lepton asymmetry which can be converted through electroweak effects to a baryon 
asymmetry, a process known as leptogenesis 80 . Furthermore the existence of an inter¬ 


mediate scale makes gauge coupling unification feasible in the absence of supersymmetry. 

Here we studied in detail one of the often unheralded features of SO(IO) grand uni- 
hcation. Namely its ability to provide for a WIMP dark matter candidate in addition 
to the benefits described above. Since SO(IO) includes an additional U(l) symmetry be¬ 
yond hypercharge, if the Higgs held that breaks this symmetry at the intermediate scale 
belongs to a 126 dimensional representation, then a discrete Z 2 symmetry is preserved 
at low energies. This discrete symmetry (equivalent to matter parity Pm = (—1)^^^“^^) 
naturally provides the stability for dark matter candidates. We considered all possible 
intermediate gauge group with broken SU(5). 

Stable SO(IO) scalar (fermion) DM candidates must be odd (even) under the Z 2 sym¬ 
metry. Therefore fermions must originate in either a 10 , 45 , 54 , 120 , 126 , 210 or 210 ' 
representation, while scalars are restricted to either a 16 or 144 of SO(IO). These multi- 
plets must be split and we gave explicit examples of hne-tuning mechanisms in order to 
retain a 1 TeV WIMP candidate which may be a SU(2 )l, singlet, doublet, triplet, or quar¬ 
tet with or without hypercharge. Fermions which are SU(2 )l singlets with no hypercharge 
are not good WIMP candidates but are NETDM candidates and these were considered 
elsewhere 24 . Our criteria for a viable dark matter model required: gauge coupling unifi¬ 


cation at a sufficiently high scale to ensure proton stability compatible with experiment; a 
unification scale greater than the intermediate scale; and elastic cross sections compatible 
with direct detection experiments. The latter criterion often requires additional Higgs 
representations to split the degeneracy of the fermionic intermediate scale representations 
if DM is hypercharged. 

Despite the potential very long list of candidates (when one combines the possible 
different SO(IO) representations and intermediate gauge groups), we found only a handful 





of models which satished all constraints. Among the scalar candidates, the F = 0 singlet 
and V = 1/2 doublet (often referred to as an inert Higgs doublet (^) are possible 
candidates for SU(4)c (8 SU(2 )l ® SU(2)r and SU(3)c SU(2)l ® SU(2 )h ® U(1)b-l 
( with or without a left-right symmetry) intermediate gauge groups. These originate from 
either the 16 or 144 of SO(IO). The latter group (without the left-right symmetry) is 
also consistent with a state originating from the 144 being a triplet under SU(2)/j. To 
avoid immediate exclusion from direct detection experiments, a mass splitting of order 
100 keV implies that the intermediate scale must be larger than about 3 x 10“® Mqut for 
a nominal 1 TeV hypercharged scalar DM particle. Some of these models imply proton 
lifetimes short enough to be testable in on-going and future proton decay experiments. 

The fermion candidates were even more restrictive. Models with V = 0 must come 
from a SU(2)2, triplet (singlets are not WIMPs). In this case only one model was found 
using the SU(4)c ® SU(2 )l ® SU(2 )i^ intermediate gauge group and requiring additional 
Higgses (already present in i?i) at the intermediate scale. Models with H = 1/2 doublets 
were found for 811(4)^; ® SU(2)i 0 U(l)/j with a singlet fermion required for mixing, and 
SU(4)(7 (g) SU(2)2, (g) SU(2)/^ with a triplet fermion for mixing. In both cases, additional 
Higgses from are required at the intermediate scale. More possibilities can be found 
if the additional Higgs are taken outside Ri. 

SO(10) almost always involves rather large representations (at least when compared 
with minimal SU(5) for example). SO(IO) itself can be broken by either a 45, 54 or 210 
representation (i?i) which determines the intermediate scale gauge group and a 126 is 
needed to break down to the SM (and preserve the needed Z 2 symmetry). The SM Higgs 
originates from a 10 and matter helds reside in three copies of 16’s. One additional rep¬ 
resentation is needed to account for DM. We have delineated the possible representations 
and necessary intermediate gauge groups needed to account for WIMP-like dark matter, 
proton stability, and gauge coupling unihcation. Some of these models are accessible for 
experimental tests. 
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Appendix 

A Dark matter candidates in SO(10) 

Here we give a group theoretical argument to classify possible DM candidates in SO(IO) 


models. We basically follow the notation of Ref. 36 in the following discussion. See also 
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Refs. 1^1^. 

Since SO(IO) is a rank-five group, we have five linearly independent Cartan generators. 
We denote them hy Hi {i = 1,... ,5). In the dual basis, they are expressed in terms of 
five-dimensional vectors as follows: 


«. = jll 2 2 1 1] , 

//s = i|0 0 1 1 1] , 

//4= h-2 0 3 - 1 1]. 

D 

ifs = [2 0 2 1 - 1] . 


(16) 


Here, Hi and H 2 are the SU(3)c Cartan generators. H 3 and H 4 are the weak isospin and 
hypercharge, T^l and Y, respectively. H 5 is given by H^ = — L) + AY. 

Every component of an SO(10) multiplet is specihed by a weight vector which is 
expressed by a set of hve integers called Dynkin labels as n = (pi,..., Jl^). Its eigenvalues 
of Hi are given by 

5 

^ ^ hjjfij , (17) 

with Hi = [hii,..., hi^]. 

The DM particle should have zero eigenvalues of Hi, H 2 , and Q = H 3 + H^. This 
condition is satisfied by the following set of weight vectors characterized by two integers 
N and M: 

= {-N N - M - N + M M) . (18) 

The hypercharge and B — L charge of the weight vector are 


Y{iin,m) = ^{N-M) , 


B — — N 


(19) 


We hnd that the N ^ M cases correspond to the hypercharged DM candidates. For 
N = M, the weight vector agrees to hn discussed in Ref. |^ . 

A convenient way to determine the matter parity of the DM candidates is to 

consider congruence classes of SO(IO) (^. Two representations R and R' of SO(IO) 
belong to the same congruence class if the difference between any weights of R and R' is 
always written as a linear combination of simple roots of SO(IO) with integer coefficients 
81 . SO(IO) has four congruence classes. Each congruence class is characterized by 
a congruence number. The congruence number of a weight vector = (pi,..., ps) is 
defined by 


c{n) = 2/ri 2/^3 - /r 4 /i 5 (mod. 4) 


( 20 ) 
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For have 


c{ij^n,m) = 3iV + 2M (mod. 4) . 


( 21 ) 


Then, we hnd that the congruence number of iin,m is ±1 if and only if N is an odd 
integer. Considering Eq. (19), we conclude that a DM particle has the odd matter parity 
if it has a congruence number of ±1, while it has the even matter parity if its congruence 
number is 0 or 2. 


B Proton decay calculation 


In this section, we describe how we calculate proton decay lifetimes in the intermediate- 
scale scenarios. In these scenarios, proton decay is induced by the exchange of the GUT- 
scale gauge bosons 82 . The relevant part of the SO(IO) gauge interaction is given by 


C-mt — 


5'gut r/7T 
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[{Q)arr"PR{L% + + h-C-] , 

( 22 ) 


where 



X represents the GUT gauge bosons which induce proton decay, (^gut is the unihed 
gauge coupling constant, a, h, c are SU(3)c indices, i,j are SU(2 )l indices, r, s are SU(2)ij 
indices, and Pr/l = (1 ± 75)/2 are the chirality projection operators. The exchange 
of the X helds generates dimension-six proton decay operators. These operators are 
expressed in a form that respects the intermediate gauge symmetries. Between the GUT 
and intermediate scales, the renormalization factors for the effective operators are in 
general different among the choices of Gint- Below the intermediate scale, the low-energy 
effective theory is described by the SU(3)c®SU(2)i®U(l)y gauge theory, and thus after 
matching the theories above and below the intermediate scale, the prescription for the 
calculation is common to all of the cases. For this reason, we hrst describe the calculation 
below the intermediate scale. After that, we discuss each intermediate gauge theory 
showing the matching conditions at the GUT and intermediate scales as well as the RGBs 
between them. 

In the SU(3)c ® SU(2)i (g) U(l)y gauge theory, the effective Lagrangian for proton 
decay is generically written as 


>Ceff = 5^G,G, + h.c. , (24) 

/=i 
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with the effective operators given by [83|-[8h 


Oi — eabcQj(ujid^)(Q'^L-l) , 

O3 = eabce^MQLQ''L){Q^L{) , 

C>4 = ea6c(MR<^R)(M^e^) , (25) 


up to dimension six. We then run down the coefficients to the electroweak scale. We will 
see below that the coefficients and vanish in all of the cases we consider in this 
paper, and thus we focus on Ci and C 2 - Their renormalization factors are 
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CM 


a3(p) 

^3 

a2(p) 

4^2 

ai(p) 

0(3 (Mint). 
03 (P) 

_ ^ 
^3 

.«2(Tfint). 

a2(p) 

_ ^ 

4^2 

_ai(Mint). 

ai(p) 

as (Mint). 


.02 (Mint). 


_Q!l(Mint). 


11 

2 O 61 


23 

20bi 


Cl(Afint) , 

C2(AAnt) , 


( 26 ) 

( 27 ) 


where ba denote the one-loop beta-function coefficients for the gauge couplings ga and /r is 
an arbitrary scale. We need to change the beta-function coefficients appropriately when 
we across the DM mass threshold. Below the electroweak scale, the QCD corrections are 
the dominant contribution. By using the two-loop RGE given in Ref. [^, we compute 
the Wilson coefficients at the hadronic scale /Xhad as 


Ci (p-had) 


(/^had) 

6 

25 

as{mb) 

6 

23 

Crs(Phad) + 

173 

825 

' asimb) + ^ ' 

_ asirub) . 


_as{mz)_ 


. as{mb) + ^ . 


_as{mz) + 


430 

2001 


Ci{mz) 




( 28 ) 


with i = 1,2. 

In non-SUSY GUTs, the dominant decay mode of proton is p —)■ Tr^e’*'. The partial 
decay width of the mode is computed as 


r(p —)■ TT^e'*') 



( 29 ) 


where rup and are the masses of the proton and the neutral pion, respectively, and 

Al = Gi(phad)( 7 r° I Ip) , 

= 2G2(phad)(7r°|(Md)LMij|p) . (30) 

The hadron matrix elements are evaluated with the lattice QGD simulations in Ref. (^ . 
We have 


{'K^\{ud)RUL\p) = {'K^\{ud)LUR\p) = -0.103(23)(34) GeV^ , (31) 

with phad = 2 GeV. Here, the first and second parentheses indicate statistical and sys¬ 
tematic errors, respectively. 
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B.l Gint = SU(4)c ® SU(2 )l 0 SU(2)fl(0D) 

For Gint = SU(4)c ® SU( 2 )l 0 SU(2)ij(0D), the dimension-six effective operator is given 

b}H _ _ 

+ h-c. , (33) 

where a,f3,... denote the SU(4) indices, and the Dirac field 4/ = is defined by 


= 


ul ul ul Vl 

,4 4 4 eL, 




4i 


42 
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'R 


'^Rl '^R2 '^R'i ^R; 


(34) 


Here, the indices represent the SU(3)c color and C indicates charge conjugation. At tree 
level, the coefficient of the effective operator is evaluated a^ 


G422(AfGUT) — 


4ut 

'2M2. ’ 


(35) 


with Mx the mass of the heavy gauge field X. In this paper, we neglect fermion flavor 
mixings for simplicity. 

The Wilson coefficient is evolved down to the intermediate scale using the RGE. The 
renormalization factor is computed to be 


G422(Tfint) — 


0^4 ( A/jnt) 

15 

464 

tt 2 L(Tfint) 

9 

tt 2 ij(Tfint) 

CKgUT 


(UGUT 


cigut 


G422(TfGUT) 


(36) 


Then, the effective operator is matched onto the operators in Eq. (|25|). The Wilson 
coefficients Cj are given b} 


16 


Gi(Mint) = 4G422(Mi„t) , 

G2(Mint) = 2G422(Mint) , 

G3(Mint) = G4(Mi„t) = 0 . (37) 


B.2 Gint = SU(4)c 0 SV{2)l 0 V{1)r 

In the case of Gint = SU(4 )g 0 SU( 2 )£, 0 U(l)ii;, the effective Lagrangian is written as 


Ces = G 421 ■ 2 e,,eo^f,^si^<^rPL^^^iU^yPRV' + h.c. , 

Note that 




and thus the operator in Eq. ( |33[ ) is the unique choice. 

^®We have found that the sign of this equation is opposite to that 


16 


given in 


Ref. 


We have fixed an error in the matching conditions given in Ref. 


24 


24 


(38) 

( 32 ) 
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with 


U = (m^, v) , V = {d^, d^, e) . 

The GUT-scale matching condition for the operator is 


G42i(Mgut) — 


Q;4(A/int) 

15 

454 

« 2 L(Tfint) 

9 


cxgut 


ckgut 


ckgut 


3 


9gvt 
'2Ml ’ 

and the renormalization factor is given by 

G42l(Mint) = 

For the intermediate-scale matching conditions, we have 

= 4^421 , 

G 2 (Mint) = 2G42l(Mint) , 

G3(Mint) = G 4 (Mi„t) = 0 . 

B .3 Gi„, = SU( 3 )c ® SU( 2 )i ® SU( 2 )b ® U(l)B_t(®n) 


(39) 


(40) 


G42i(Mgut) • (41) 


(42) 


When Gint = SU(3 )g G) SU( 2 )i (g) SU( 2 )/j (g) U(1 )b-l(®-D), there are four independent 
effective operators (^ , 

Qi = 2ei,erseabc{WrPLQ^KWY''PRL^ , 

Q2 = 2ei,erseabYWrPLLYWtPRQ^^ , 

Qs = 2eae,keabc{WrPLQ^KWr''PLL^ , 

Q4 = 2epse,reabc{WrPRQ"YWY^PRL^ , 

and thus the effective Lagrangian is expressed as 

4 

Geflf = ^ G3221 Ql + h.C. 


(43) 


(44) 


7=1 


For the GUT-scale matching condition, we have 


C^^22li.^GVT) - C32L(^GUt) - 


C'32L(^GUt) — G322i(Mgut) — 0 . 

The renormalization factors for the coefficients G 3221 and G 3221 are given in Refs. 


(45) 
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90 


G(Mi,t) 

G3(^int) 

2 

^3 

«2L(Tfint) 

9 

462T 

tt2R(Tfint) 

9 

aB-I,(Tfint) 

C(Mgut) 

ggut 


ckgut 


ggut 


(UguT 


46i; 


(46) 
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forC = cSi andCgi. 


Then the Wilson coefficients at the electroweak scale are matched 


onto those of the operators (25) as 


Cl (Mint) = 4C'gi(Mint) , 

C'2(Mint) = 2C'g)i(Mint) , 

C3(Mint) = C' 4 (Mi„t) = 0 . (47) 


C Example of fine-tuning 


To show the process of mass hne-tuning explicitly, in this section, we consider the case of 
-Rdm = 16 with Gint = SU(3)c®SU(2)L(g)SU(2)ii;®U(l)B_L as an example. We take = 
45, which contains two independent SM singlet components that might develop VEVs; one 
is in a (1,1, 3) while the other is in a (15,1,1) under SU(4)(^(8)SU(2)^(8)SU(2)j:j. We refer 
to these VEVs as Ai and ^ 42 , respectively, and other notation is taken from Eq. ([^. Since 
the components of a scalar 16 have the same quantum numbers as those of a generation 
of the SM fermions, we denote them by the same symbol as for the corresponding SM 
fermions with a tilde, just like the notation for sfermions in supersymmetric models. 

Let us hrst study the coupling. Since Ri is the adjoint representation of 

SO(IO), the decomposition of this coupling in terms of the component helds has a similar 
form to the gauge interaction for a 16 spinor representation. We have 




— \/3y42^ T (^y^Ai — \/3A2^ Rr^r T '/^A2R*^L 
( H—^^2 ] (RrcIj^ + ( — y/2Ai H—^^42 J - -^A2Q*j^Q^ 
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73 


73 


(48) 


where the contraction of the SU(3)c and SU( 2 )i indices is implicit. When 7 ^ 0 and 
A 2 = 0, the mass spectrum preserves the SU(4)(^ 0 SU(2)^ 0 U(l)^ symmetry, while 
when ^2 7 ^ 0 and Ai = 0, then it is SU(3)(^ 0 SU(2)^ 0 SU(2)^ 0 U(l)^_^ symmetric. If 
both of the VEVs have non-zero values, then the low-energy theory is invariant under the 
SU(3)(^ 0 SU(2)^ 0 U(l)^ 0 U(l)^_^ symmetry. The coefficients of A 2 for left and right 
doublets have different signs, which indicates the breaking of left-right symmetry. Here, 
we choose Hi = 0 and H 2 = V 45 to obtain Gint = SU(3)c 0 SU( 2 )i 0 SU( 2 )ii; 0 U( 1 )s_l. 

Next we consider the mass terms generated by Af^(i?oMi?D]y[) 45 (i? 2 -R 2 ) 45 - The SM 
singlet in R 2 = 126 transforms as (10,1,3) under SU(4)(^ 0 SU(2)j- 0 SU(2)j:j, which 
acquires a VEV ni 26 to break Gint into the SM gauge group. According to the result in 
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Ref. 91 , the resultant mass terms an 


Af (flDMflDM)45 mR2U) = 


^126 


'^R^R + g ( + d*Rdj^ 


g ( ^R^R + '^r'^R + QlQl 


(49) 


Notice that the right-hand side of the expression can be grouped in terms of SU(5) multi- 
plets. This is expected since is invariant under the SU(5) transformation. From the 
above equations, it is found that we can ensure that only the DM component has a mass 
around TeV scale by fine-tuning the parameters ki and Af^. For example, to obtain 
the model SA 3221 , we can take 


m 2 - y3«in45 ~ 0{Ml,) , 

m 2 - y3Kin45 - Af ~ 0(TeV2) . (50) 

Then, Dr acquires a TeV-scale mass, while the mass of Dr is The rest of the 

components lie around the GUT scale. On the other hand, if we take 

m 2 + V3«:in45 ~ 0{Ml,) , 

m2 + ^ 3 ^ 1^45 + ?Af ^2^6 ~ C>(TeV2) , (51) 

5 

then we can make only the Lr component have a TeV-scale mass and the other components 
have GUT-scale masses. Thus we obtain the SB 3221 model. 

To simplify our argument, in the above discussion, we have taken into account only 
the contribution of the M2, ki, and Af terms, and neglected that of the other terms in 
Eq. ([^. Even in the presence of the other contributions, we can always perform a similar 
hne-tuning among the parameters to realize desired mass spectrum for our DM models. 


^^Note that since (i? 2 -^ 2)45 contains a 45, there is a contribution to the mass corresponding to Eq. (48) 
at the intermediate scale proportional to Af® with independent coefficients Ai and 442. The result shown 
is obtained from Eq. (48) by taking Ai = and 442 = up to an overall factor. 
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